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Abstract 

Background: The reproductive phenology of perennial plants in temperate climates is largely conditioned by the 
duration of bud dormancy, and fruit developmental processes. Bud dormancy release and bud break depends on 
the perception of cumulative chilling and heat during the bud development. The objective of this work was to 
identify new quantitative trait loci (QTLs) associated to temperature requirements for bud dormancy release and 
flowering and to fruit harvest date, in a segregating population of peach. 

Results: We have identified QTLs for nine traits related to bud dormancy, flowering and fruit harvest in an 
intraspecific hybrid population of peach in two locations differing in chilling time accumulation. QTLs were located 
in a genetic linkage map of peach based on single nucleotide polymorphism (SNP) markers for eight linkage 
groups (LGs) of the peach genome sequence. QTLs for chilling requirements for dormancy release and blooming 
clustered in seven different genomic regions that partially coincided with loci identified in previous works. The 
most significant QTL for chilling requirements mapped to LGl, close to the evergrowing locus. QTLs for heat 
requirement related traits were distributed in nine genomic regions, four of them co-localizing with QTLs for chilling 
requirement trait. Two major loci in LG4 and LG6 determined fruit harvest time. 

Conclusions: We identified QTLs associated to nine traits related to the reproductive phenology in peach. A search 
of candidate genes for these QTLs rendered different genes related to flowering regulation, chromatin modification 
and hormone signalling. A better understanding of the genetic factors affecting crop phenology might help 
scientists and breeders to predict changes in genotype performance in a context of global climate change. 
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Background 

The timing of reproductive events of perennial plants in 
temperate climates is largely conditioned by dormancy, a 
period of cyclic quiescence during the low temperatures 
of autumn and winter, which protects meristems within 
buds from the detrimental effects of cold and water 
stress. The pioneering work by Lang [1] distinguished 
the dormancy due to mechanisms intrinsic to the bud 
(endodormancy) from the inhibition of growth imposed 
by other organs of the plant (paradormancy) or by envir- 
onmental factors (ecodormancy). However, more recent 
reviews tend to emphasize the idea of dormancy as a 
state within the meristem, independently of the origin of 
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the dormancy-imposing cues; and highlight the dynamic 
and quantitative nature of dormancy, varying in intensity 
according to intrinsic and environmental signals [2,3]. 

The quantitative perception of environmental chilling 
is the major and best-known factor favouring the release 
of bud dormancy [4]. The extent of chilling needed is 
highly genotype-dependent and constitutes an adaptive 
strategy to the duration of the cold season under spe- 
cific climate conditions. After fulfilment of the chilling 
requirements for dormancy release, a period of warm 
temperatures is needed prior to bud burst (heat require- 
ment) and the subsequent developmental phases leading 
to fruit set, growth and maturation. Both, chilling re- 
quirement for dormancy release and fruit maturation 
time have been described as two major limiting factors 
determining respectively the southern and northern 
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boundaries to the geographical distribution of temper- 
ate species [5]. 

Recently researchers have increased their interest on 
the effects of global climate warming on plant phenology 
[6]. Whereas some species are expected to advance their 
growing season due to the increasing temperatures in 
spring, others could delay or experience abnormal bud 
burst as a result of an insufficient chilling for dormancy 
release in winter. In particular, a delay in the beginning 
of the growing season observed in the meadow and 
steppe vegetation of the Tibetan Plateau from the mid- 
1990s has been related to the later fulfilment of chilling 
requirements [7]. 

Plant species cope with changing climate conditions 
by shifting their geographical distribution, with a plastic 
response of plant phenology to environmental changes, 
or alternatively through the natural selection of popula- 
tions with dormancy, flowering and fruit maturation 
traits adapted to the new conditions [8]. The knowledge 
on the genetic factors affecting these phenology-related 
traits is scarce and fragmentary in woody perennial spe- 
cies, although recent remarkable advances have been 
reached based on the comparison with analogous pro- 
cesses of annual model species, the use of novel 
transcriptomic and molecular approaches, and the ex- 
ploitation of the natural variability by means of mutant 
and quantitative trait locus (QTL) analyses. Some of the 
most relevant QTL studies on reproductive phenology 
have been conducted in poplar and in species within the 
Rosaceae family. In poplar these studies have focused 
mostly on bud set and bud flush in spring, resulting in 
some loci co-localizing with genes involved in light per- 
ception and abscisic acid (ABA) signalling [9-11]. In 
apple, QTL analysis of bud break date in flower and 
vegetative buds of two different progenies pointed to 
two major genomic regions controlling these traits in 
linkage group (LG) 8 and LG9, showing numerous genes 
involved in cell cycle control [12]. Blooming date and 
fruit maturation date traits have been frequently ana- 
lyzed in QTL analyses of species from the genus Prunus, 
such as almond, peach, apricot and sweet cherry [13-21]. 
In order to better characterize at the genetic level the 
physiological response of buds to chilling, some QTL 
studies included the chilling and heat requirement traits 
in their analyses [22-24]. In spite of the numerous genes 
proposed as putative candidates for the different QTLs de- 
scribed in these reports, none of them has been function- 
ally involved in dormancy or flowering regulation, with 
the exception of DORMANCY-ASSOCIATED MADS- 
BOX {DAM) genes [25,26]. 

The objective of this work was to identify new QTLs 
and candidate genes associated to chilling and heat re- 
quirements in a Prunus persica [L.] Batsch (peach) 
population segregating from a sequential cross of three 



varieties with different dormancy behaviour and origin. 
The fruit harvest time trait was also analyzed in order to 
gain insight into the effect of dormancy duration on the 
subsequent fruit phenology. 

Results 

Phenotypic assessment and correlation between traits 

Traits evaluated in the 'V6' x 'Granada' progeny are 
listed in Table 1. The mean chilling requirement values 
for dormancy release, measured following the most 
popular Weinberger (CRW), Utah (CRU) and dynamic 
(CRD) models, were essentially equivalent in AA and EJ, 
in spite of the two-weeks delay of chilling accumulation 
observed in EJ location with respect to AA at 500 CU 
and 500 CH (Figure 1). This suggests that the three 
models used for the evaluation of chilling requirements 
are essentially valid under the climatic conditions of AA 
and EJ sites. In close agreement with the chilling delay 
of two weeks in EJ, the mean endodormancy release date 
in this location was 15 days later than in AA (Table 1). 
However the mean ecodormancy release and blooming 
delays between EJ and AA decreased to 11 and 8 days 
respectively, due to the faster fulfilment of heat require- 
ments in the warmer conditions of EJ location. 
Favourable temperatures also accounted for the earlier 
fruit harvest in EJ with respect to AA. On average, the 
time from dormancy release to flower and fruit develop- 
ment was 23 days shorter in EJ than in AA. 

Variables CRW-AA, CRW-EJ, CRU-EJ, CRD-EJ, HRB-AA, 
EcD-AA, BD-AA, PEnB-AA and HD-2012-AA departed 
from normality due to altered standardized skewness 
or kurtosis of their frequency distribution (Figure 2). 
CRW-AA and CRW-EJ distributions were particularly 
skewed to the left, that is enriched in low chilling 
individuals. EcD and PEnB distributions showed bimodal 
or multimodal profiles, with two or more separated peaks 
in both locations. 

The parental cultivars of 'V6' selection, 'Maruja' and 
'Red Candem', had the highest and lowest range values 
of chilling requirement distributions (CRW, CRU and 
CRD) and date variables (EnD, EcD, BD and HD) re- 
spectively. 'Maruja' and 'Red Candem' data were not 
available for variables HD-2011-EJ and HD-2012-AA. 
However traits related to heat/ time requirements for 
ecodormancy release and blooming (HREc, HRB, PEnEc 
and PEnB) showed a more variable location of parental 
genotypes within the range, and numerous transgressive 
segregants exceeded the parent values (Figure 2). 

The Pearson correlation coefficients between variables 
are shown in Table 2. The chilling requirement (CRW, 
CRU and CRD), EnD, EcD and BD traits showed strong 
positive correlations among them (r>0.84, P<0.01), 
moderate positive correlations with HD (r > 0.30, P < 
0.05), and strong negative correlation with PEnB. Heat 



Romeu ef al. BMC Plant Biology 2014, 14:52 
http://www.biomedcentral.com/1471-2229/14/52 



Page 3 of 16 



Table 1 Traits investigated in the 'V6' x 'Granada' progeny, with the mean value, standard deviation and data range 
over the whole population in locations AA and EJ 



Trait 


unit 


ADD. 


AA population 


EJ population 


M6an 


SD 


R3nc|6 


M6an 


SD 




ChillinQ recjuirement for GndodormBncy tgIgssg 


















Weinberger model 


Chilling hours (CH) 


CRW 


400 


118 


248-760 


373 


138 


212-658 


Utah model 


Chilling units (CU) 


CRU 


577 


154 


387-960 


591 


145 


380-856 


Dynamic model 


Portions 


CRD 


37 


8 


25-58 


42 


8 


30-58 


Endodormancy release date 


Julian days 


EnD 


20 


12 


1-51 


35 


14 


15-66 


Ecodormancy release date 


Julian days 


EcD 


46 


13 


23-64 


57 


12 


36-73 


Blooming date 


Julian days 


BD 


67 


8 


58-82 


75 


8 


57-93 


Heat requirement for ecodormancy release 


growing degree hours 
(GDH) 


HREc 


2500 


627 


908-4004 


2797 


835 


1 1 05-4670 


Heat requirement for blooming 


growing degree hours 
(GDH) 


HRB 


5110 


669 


2801-7372 


5956 


830 


3826-7956 


Period of time between endo- and ecodormancy release 


Days 


PEnEc 


26 


6 


11-41 


22 


7 


7-37 


Period of time between endodormancy release and 
blooming 


Days 


PEnB 


46 


7 


25-60 


40 


8 


23-58 


Harvest date 


Julian days 


HD 


187 


19 


155-218 


171 (2011) 
179 (2012) 


13 
14 


144-192 
1 50-204 



Abb., variable abbreviation; SD, standard deviation. 

requirements traits (HREc and HRB) were also correlat- 
ing with their respective time interval counterparts 
(PEnEc and PEnB) with r > 0.60 and P < 0.01. Correla- 
tions between AA and EJ locations were high for most 
of the traits, with the exception of HREc (r = 0.15, P = 
0.24), HRB (r=0.23, P = 0.06) and PEnEc (r = 0.24, P = 
0.05) (Table 3), which indicates a strong genotype-x- 
environment interaction in these latter traits. 



Map construction 

Genotyping was performed using the International 
Peach SNP Consortium (IPSC) peach 9 K Infinium" II 
array [27]. Briefly, 2,865 SNPs from the array were iden- 
tified as polymorphic (40% of the total) and used for 
map construction. Since the different linkage groups 
showed several SNPs co-localizing at the same locus, 
one SNP per locus was selected in order to obtain a 




-40 -20 0 20 40 60 

Date (Julian days) 



Figure 1 Chilling accumulation in AA and EJ locations along the cold season. Chilling accumulation in AA (circles) and EJ (triangles) 
locations, according to Weinberger (white symbols) and Utah (black symbols) models. The discontinuous line is to interpolate the date with a 
chilling accumulation value of 500 CH/CU. 
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Figure 2 Frequency distribution of traits measured in the 'V6' x 'Granada' progeny. Chilling requirement (CRW, CRU and CRD) and heat 
requirement (HREc and HRB) variables are represented in A. The values for the V6' (V6), 'Granada' (G), 'Red Candem' (RC) and 'Maruja' (M) parents 
are indicated by arrows. Endodormancy release date (EnD; black bars), ecodormancy release date (EcD; grey bars), and blooming date (BD; white 
bars) are represented in B. The arrows indicating date values of parents are black for EnD, grey for EcD and white for BD. The panel C includes 
the endodormancy-ecodormancy (PEnEc) and endodormancy-blooming (PEnB) time periods. Harvest date (HD) trait is represented in panel D. 



condensed map. For the 'V6' map, a total of 178 SNPs 
were retained that covered all the chromosomes but 
chromosome 2, representing a total distance of 480 cM 
with an average marker density between adjacent 
markers of 2.94 cM. For 'Granada' map, 76 SNPs were 
retained, covering chromosomes 2, 4, 5, 6, 7 and 8 with 
a total distance of 276 cM and an average marker dens- 
ity of 3.87 cM. This lack of polymorphic markers in 
large chromosome regions should be due to homozygos- 
ity of those regions. Further details on genetic map con- 
struction can be found in [28]. 



QTL analysis of chilling requirement and flowering time 

The analysis of co-segregation between SNPs and 
phenotypic data led to the identification of QTLs for all 
the investigated traits. For every variable, the LOD 
threshold for P < 0.05 calculated by the permutation test 
was between 2.2 and 2.5. Several major QTLs explaining 
60-76% of the phenotypic variance of CRW, CRU, CRD, 
EnD, EcD and BD overlapped within the genomic region 
lb, close to SNP_IGA_1 22057 marker at the end of LGl 
in 'V6' map (Table 4, Figure 3). These QTLs in lb were 
consistently detected in the AA and EJ locations. 
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Table 2 Pearson correlation coefficients between variables in the 'V6' x 'Granada' progeny 



Variable Location CRW CRU CRD EnD EcD BD HREc HRB PEnEc PEnB HD-2011 



CRU 


AA 


0.99" 






















EJ 


0.98** 




















CRD 


AA 


0.99** 


1 .00" 




















EJ 


0.98** 


1 .00" 


















EnD 


AA 


0.99" 


1 .00** 


1 .00** 


















□ 


0.98" 


1 .00** 


1 .00** 
















EcD 


AA 


0.87" 


0.90** 


0.90** 


0.90** 
















EJ 


0.84" 


0.88" 


0.88** 


0.87** 














BD 


AA 


0.86" 


0.86" 


0.86** 


0.86** 


0.85** 














EJ 


0.87" 


0.87" 


0.86** 


0.87** 


0.88** 












HREc 


AA 


0.31" 


0.32" 


0.32** 


0.31** 


0.59** 


0.43** 












EJ 


0.11 


0.12 


0.11 


0.09 


0.52** 


0.36** 










HRB 


AA 


-0.01 


-0.01 


-0.02 


-0.03 


0.08 


0.46** 


0.33** 










EJ 


-0.09 


-0.13 


-0.14 


-0.14 


0.11 


0.36** 


0.58** 








PEnEc 


AA 


-0.03 


0.03 


0.03 


0.02 


0.45** 


0.19 


0.74** 


0.25* 








EJ 


-0.55" 


-0.50" 


-0.50** 


-0.52** 


-0.03 


-0.23* 


0.72** 


047** 






PEnB 


AA 


-0.80** 


-0.80" 


-0.80** 


-0.81** 


-0.64** 


-0.39** 


-0.05 


0.60** 


0.18 






EJ 


-0.83" 


-0.85" 


-0.85** 


-0.85** 


-0.61** 


-0.48** 


0.22* 


0.63** 


0.68** 




HD-2011 


EJ 


0.37" 


0.36** 


0.37** 


0.36** 


0.36** 


0.40** 


0.22 


0.14 


-0.07 


-0.22 


HD-2012 


AA 


0.30* 


0.30* 


0.31* 


0.32** 


0.44** 


0.32** 


0.31** 


-0.01 


0.36** 


-0.21 




EJ 


0.37" 


0.38** 


0.39** 


0.37** 


0.40** 


0.37** 


0.24* 


0.06 


-0.04 


-0.26* 0.87** 



*P < 0.05, **P < 0.01. 



Another cluster of QTLs was identified in the LG7 of 
'Granada' map. The LOD plots showed two close linked 
peaks, suggesting the possibility that two linked QTLs 
were present in this genomic region (Table 4). Given 
that the 2-LOD confidence intervals of these putative 
linked QTLs were shortly overlapping (Table 4) we con- 
sidered that region as a single QTL, although the possi- 
bility of two linked QTLs could not be discarded. This 



Table 3 Pearson correlation coefficients between 
variables measured in different locations (AA/EJ) 



Variable 


AA/EJ 


P-value 


CRW 


0.88 


<0.01 


CRU 


0.89 


<0.01 


CRD 


0.89 


<0.01 


EnD 


0.89 


<0.01 


EcD 


0.87 


<0.01 


BD 


0.79 


<0.01 


HREc 


0.15 


0.24 


HRB 


0.23 


0.06 


PEnEc 


0.24 


0.05 


PEnB 


0.61 


<0.01 


HD-2012 


0.96 


<0.01 



region in LG7 contained QTLs for the traits chilling re- 
quirement, EnD and BD, which were consistently de- 
tected in the two locations with ranging between 14% 
and 29%. Also a QTL for EcD was detected in the AA 
location explaining 30% of the genetic variance. Other 
QTLs with minor effects (r < 18%) clustered in LG3 
and LG6, although they were detected in only one loca- 
tion (Table 4, Figure 3). 

QTL analysis of heat requirement and ecodormancy 
release 

Traits describing processes following endodormancy re- 
lease were HREc and HRB, which use the Anderson 
model to estimate the heat requirements for ecodor- 
mancy release and blooming respectively. Two add- 
itional traits, PEnEc and PEnB, estimate the time in days 
required for the fulfilment of HREc and HRB values re- 
spectively, starting from the endodormancy release date. 
A QTL cluster for HREc and PEnB traits with consistent 
major effects in both locations (R = 15-39%) was de- 
tected in LGl of 'V6' map, overlapping with QTLs for 
chilling requirement and flowering time described above. 
Additional QTLs with minor effects but not consistent 
across locations were also detected in LG3, LG4, LG5, 
LG6 and LG7 (Table 4, Figure 3). 
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Table 4 Quantitative trait loci (QTLs) detected in the 'V6' x 'Granada' progeny 



Variable 


Location/year 


Map 


LG 


Position (cM) 


Nearest marker 


CI (cM) 


LOD 


(%) 


Additive effect 


CRW 


AA 


V6 


1 


72.5 


SNP_IGA_1 22057 


68.2-74.5 


16.3 


69 


-206.6 




EJ 


V6 


1 


73.5 


SNP_IGA_1 22057 


68.2-74.5 


22.3 


64 


-228.8 




EJ 


V6 


3 


104 


SNP_IGA_297497 


8.9-14.4 


3.8 


10 


-123.2 




EJ 


G 


7 


29.9 


SNP_IGA_769194 


29.3-40.5 


3.9 


20 


-159.7 




AA 


G 


7 


45.5 


SNP_IGA_779224 


41.3-48.9 


3.9 


24 


-120.1 


CRU 


AA 


V6 


1 


72.5 


SNP_IGA_1 22057 


66.8-74.5 


16.3 


76 


-284.2 




EJ 


V6 


1 


73.5 


SNP_IGA_1 22057 


69.3-74.5 


22.1 


72 


-248.1 




EJ 


G 


6 


34.2 


SNP_IGA_635355 


23.6-37.1 


3.8 


18 


-143.8 




EJ 


G 


7 


29.9 


SNP_IGA_769194 


29.3-41.3 


3.7 


18 


-158.5 




EJ 


G 


7 


45.5 


SNP_IGA_779224 


43.2-48.9 


4.5 


21 


-151.2 




AA 


G 


7 


45.5 


SNP_IGA_779224 


41.3-48.9 


3.9 


25 


-158.5 


CRD 


AA 


V6 


1 


72.5 


SNP_IGA_1 22057 


66.8-74.5 


17.3 


65 


-141 




EJ 


V6 


1 


73.5 


SNP_IGA_1 22057 


68.4-74.5 


22.1 


67 


-13.9 




EJ 


V6 


3 


2.9 


SNP_IGA_293752 


0.0-15.6 


2.5 


6 


-4.2 




EJ 


G 


6 


34.2 


SNP_IGA_635355 


23.3-41.6 


2.8 


13 


-7.1 




AA 


G 


7 


45.5 


SNP_IGA_779224 


41.7-48.9 


4.0 


25 


-84 




EJ 


G 


7 


45.5 


SNP_IGA_779224 


41.3-48.9 


3.2 


14 


-6.5 


EnD 


AA 


V6 


1 


72.5 


SNP_IGA_1 22057 


66.8-74.5 


17.0 


64 


-26.6 




EJ 


V6 


1 


73.5 


SNP_IGA_1 22057 


67.8-74.5 


20.2 


63 


-22.2 




EJ 


G 


6 


34.2 


SNP_IGA_635355 


23.4-41.6 


3.0 


14 


-11.9 




AA 


G 


7 


45.5 


SNP_IGA_779224 


41.9-48.9 


3.9 


25 


-124 




EJ 


G 


7 


45.5 


SNP_IGA_779224 


41.3-48.9 


3.5 


16 


-11.2 


EcD 


AA 


V6 


1 


74.5 


SNP_IGA_1 22057 


72.5-74.5 


15.9 


65 


-134 




EJ 


V6 


1 


73.5 


SNP_IGA_1 22057 


67.4-74.5 


20.2 


63 


-19.2 




AA 


V6 


4 


2.8 


SNP_IGA_381379 


0.0-12.8 


2.8 


8 


-8.1 




AA 


G 


7 


45.5 


SNP_IGA_779224 


43.3-48.9 


5.2 


30 


-15.3 


BD 


AA 


V6 


1 


72.5 


SNP_IGA_1 22057 


65.8-74.5 


16.8 


60 


-224 




EJ 


V6 


1 


72.5 


SNP_IGA_1 22057 


67.2-73.5 


19.9 


60 


-12.9 




EJ 


V6 


3 


17.8 


SNP_IGA_298293 


15.2-18.8 


5.2 


16 


-9.6 




EJ 


V6 


3 


33.1 


SNP_IGA_316315 


32.5-41.9 


3.9 


11 


8.1 




AA 


V6 


6 


43.2 


SNP_IGA_664540 


42.2-47.7 


5.1 


17 


-11.6 




AA 


G 


7 


29.9 


SNP_IGA_769194 


29.3-41.3 


2.6 


19 


-8.8 




EJ 


G 


7 


29.9 


SNP_IGA_769194 


29.3-41.3 


3.7 


22 


-9.8 




AA 


G 


7 


45.5 


SNP_IGA_779224 


43.4-48.9 


4.5 


29 


-9.5 




EJ 


G 


7 


45.5 


SNP_IGA_779224 


41.3-48.9 


3.6 


20 


-8.3 




A A 


VD 


1 
1 


/J.J 


jINr loA 1 ZZUJ / 


OJ. 1 / ^.J 


J.U 


1 ^ 
1 J 


AQ7 1 




EJ 


V6 


1 


74.5 


SNP_IGA_1 32237 


63.7-74.5 


54 


25 


-644.8 




AA 


V5 


3 


124 


SNP_IGA_297497 


0.0-16.7 


2.5 


15 


-253.7 




EJ 


G 


5 


204 


SNP_IGA_591439 


2.8-32.5 


3.0 


17 


432.3 


HRB 


AA 


V6 


6 


50.6 


SNP_IGA_579852 


49.7-54.1 


5.6 


32 


457.9 




EJ 


G 


4 


114 


SNP_IGA_513496 


0.0-20.6 


2.7 


12 


584.2 


PEnEc 


EJ 


V6 


1 


14.3 


SNP_IGA_7895 


2.8-19.3 


2.8 


15 


-5.9 




EJ 


G 


5 


154 


SNP_IGA_591439 


9.8-34.1 


2.7 


14 


-5.0 
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Table 4 Quantitative trait loci (QTLs) detected in the 'V6' x 'Granada' progeny (Continued) 



PEnB 


AA 


V6 


1 


73.5 


SNP_IGA_1 22057 


65.8-74.5 


8.8 


39 


8.6 




EJ 


V6 


1 


73.5 


SNP IGA 122057 


63.7-74.5 


8.3 


33 


9.0 




AA 


V6 


5 


9.5 


SNP IGA 556166 


0.0-16.2 


2.9 


12 


4.8 




EJ 


G 


6 


34.2 


SNP IGA 635355 


23.5-37.1 


4.3 


18 


6.5 




EJ 


G 


7 


45.5 


SNP IGA 779224 


41.3-48.9 


2.6 


9 


4.9 


HD 


AA/2012 


V6 


1 


70.8 


SNP_IGA_107581 


63.6-73.5 


7.8 


16 


-16.1 




EJ/2012 


V6 


1 


73.5 


SNP_IGA_1 22057 


65.6-74.5 


7.8 


15 


-10.9 




AA/2012 


V6 


4 


52.9 


SNP_IGA_411147 


50.5-56.9 


20.1 


51 


-27.9 




EJ/201 1 


V6 


4 


52.9 


SNP_IGA_411147 


50.5-57.0 


16.2 


49 


-19.4 




EJ/2012 


V6 


4 


53.9 


SNP_IGA_411147 


52.9-56.8 


21.1 


54 


-20.5 




AA/2012 


V6 


6 


11.1 


SNP_IGA_620767 


0.0-20.9 


2.8 


6 


10.0 




EJ/201 1 


V6 


6 


11.0 


SNP_IGA_620767 


0.0-23.9 


34 


8 


7.7 




EJ/2012 


V6 


6 


11.1 


SNP_IGA_620767 


1.0-19.9 


54 


10 


8.9 




AA/2012 


V6 


7 


4.0 


SNP_IGA_713270 


0.0-17.8 


4.0 


9 


11.7 



LG, linkage group; CI, two-LOD confidence interval of QTL position; LOD, logarithm of the odds; R^, percentage of the phenotypic variance explained by the QTL. 



QTL analysis of fruit harvest date 

QTLs for HD were detected in four genomic regions of 
LGl, LG4, LG6 and LG7 of the 'V6' map (Table 4, 
Figure 3). All of them were consistent across geographical 
locations and years, except for a minor QTL in LG7. The 
major QTL was located in region 4b, with LOD values be- 
tween 16.2 and 21.1 and around 50%. A second cluster 
of QTLs was associated to the region lb, co-localizing with 
the most significant QTL for chilling requirement (CRW, 
CRU and CRD), EnD, EcD, BD, HREc and PEnB, which 
highlights the outstanding influence of this locus on the re- 
productive phenology of these peach cultivars. Other minor 
QTLs, explaining 10% or less of the phenotypic variance of 
HD, were localized in genomic regions 6a and 7a. The most 
significant QTLs in regions lb and 4b showed negative ef- 
fects on HD. 

Search of candidate genes 

QTL regions lb, 3a, 4b, 6a, 6b and 7b, containing multiple 
QTLs and QTLs with high significance, relevance and con- 
sistent effects among trials were searched for the presence 
of plausible candidate genes. The QTL regions limited by 
the genome coordinates shown in Additional file 1 were 
examined manually. Those genes showing high similarity 
to genes involved in molecular and physiological processes 
previously associated to regulation of bud dormancy, 
flowering and fruit maturation, such as ABA regulation, 
cold acclimation, ethylene signalling, chromatin modifica- 
tion, flowering and vernalization pathways were selected as 
candidate genes (Table 5). Also the transcriptional regula- 
tors ppa008979m, ppa012329m, ppa011123m {DAM4), 
ppa010822m (DAMS) and ppa010714m {DAM6) were 
included in Table 5 given their dormancy-dependent ex- 
pression in flower buds [29-31]. These last four genes were 



added to the candidate list due to the outstanding import- 
ance of DAM genes in dormancy regulation, which were 
previously postulated as candidate genes in other QTL 
studies [22,23], and their extreme localization at the end of 
LGl after the last SNP marker flanking region lb 
(SNP_IGA_132237). 

An alternative approach for the identification of candidate 
genes has been performed by reciprocal blast analysis of 
known genes involved in chromatin modification and flow- 
ering. Reciprocal blast analysis allowed the identification of 
putative orthologs of these genes in peach located within 
the genomic regions outlined by QTL analysis. Peach puta- 
tive orthologs of genes coding for subunits of the Polycomb 
Repressive Complex 1 and 2 (PRCl, PRC2), Trithorax group 
proteins [32], histone demethylases with Jumonji (JMJ) do- 
main, histone acetyltransferases and deacetylases [33] 
(Pandey et al. 2002), and flowering factors are listed in 
Table 6. Only those putative orthologs located within a 
QTL region are shown. Some candidate genes, such as 
ppa001213m, ppa000318m, ppa000228m, ppa000162m, 
ppa005747m and ppa012369m appear in both Tables 5 
and 6; however some putative orthologs that are not located 
within one of the major QTL regions lb, 3a, 4b, 6a, 6b and 
7b, and candidate genes that have not be considered as pu- 
tative orthologs by reciprocal blast analysis appear just once. 

Discussion 

In this work we have studied nine traits related to bud 
dormancy, flowering and fruit harvest in a hybrid popu- 
lation of peach. The high level of correlation between 
traits and the clustering of QTLs in certain map posi- 
tions argue for a considerable degree of redundancy that 
recommends the joint analysis of traits in three major 
groups. Chilling requirement (CRW, CRU and CRD), 
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Figure 3 Location of quantitative trait loci (QTLs) in the 'V6' x 'Granada' genetic map. QTLs are labelled by boxes (two-LOD confidence 
interval). Boxes are black for chilling requirement and blooming variables (CRW, CRU, CRD, EnD, EcD and BD), white for heat requirement 
variables (HREc, HRB, PEnEc and PEnB) and grey for harvest date (HD). V5' and 'Granada' genetic maps are vertical white bars showing single 
nucleotide polymorphisms (SNP) markers. Genome scaffolds are represented as vertical black bars placed between their respective V6' and 
'Granada' maps. The positions of SNP markers flanking the confidence intervals of QTLs are labelled in bold on genome scaffolds, limiting a 
genomic region named with a number (for the scaffold) and a letter {for the order of the region along the scaffold), which was used to identify 
putative candidate genes. Within each genomic region, QTLs of the same variable corresponding to different locations are fused in a single box. 
Scale bars for genetic maps and scaffolds are shown in the upper left part of the figure. 
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Table 5 Peach genes located in the two-LOD confidence Interval of selected QTLs and other candidate genes 

QTL region Peach model (gene) Genomic location (Mb) Blastp hit in Arabidopsis E-value Description 



lb 


ppa000228m 


33.9 


PKL 


0 


Chromatin remodelling factor PICKLE 




ppa008143m 


34.3 


C0L4 


2x10-"° 


Zinc finger protein CONSTANS-LIKE 4 




ppa003748m 


34.8 


HABl 


8x1 0"^" 


Protein phospliatase 2C 




ppa006590m 


35.0 


HDA2 


0 


Histone deacetylase 2 




ppa008979m 


35.9 


At5g67300 


2x10* 


myb-related protein 




ppa006503m 


36.2 


AREBl 


3x10"^^^ 


ABA-responsive element binding protein 1 




ppa000318m 


37.0 


EMFl 


8x1 0-'" 


Embryonic flower 1 




ppaOl 3757m 


37.6 


FPFl 


8x10"" 


Flowering promoting factor 1 




ppa000056m 


39.8 


EPS 


0 


Histone-lysine N-methyltransferase 




ppa001566m 


40.0 


SWI3C 


0 


SWI/SNF complex subunit SWI3C 




ppa005747m 


40.8 


HAM2 


0 


Histone acetyltransferase of the MYST family 




ppa024363m 


41.0 


AREB3 


3x10-"= 


ABA-responsive element binding protein 3 




ppaOl 2329m 


46.0 


At4g24440 


3xl0-<^« 


Transcription initiation factor HA subunit 2 




ppa010714m i.DAM6) 


46.3 


SVP 


7x10"*^= 


MADS-box protein 




ppaOl 0822m [DAMS] 


46.4 


SVP 


1x10"*^ 


MADS-box protein 




ppa011123m {DAM4) 


46.4 


AGL24 


2x10"*° 


MADS-box protein 


3a 


ppa020502m 


0.1 


CCR2 


2x10-1' 


Glycine-rich protein 




ppa004975m 


0.1 


SDG40 


0 


SET domain group 40 




ppa000162m 


1.2 


HACl 


0 


Histone acetyltransferase of the CBP family 




ppaOl 3609m 


1.7 


Atlg54070 


1x10-" 


Dormancy/auxin associated family protein 




ppa004252m 


2.4 


HABl 


0 


Protein phosphatase 2C 




ppa002515m 


3.2 


AB05 


0 


Pentatricopeptide repeat-containing protein 


4b 


ppaOl 0982m 


10.4 


ERF4 


6x10-"= 


Ethylene responsive element binding factor 4 




ppa022739m 


10.5 


SPL4 


3x10-"' 


Squamosa promoter binding protein- ike 4 




ppa008301 m 


11.1 


ANAC072 


2x10-^^^ 


NAC domain-containing protein 72 




ppa002986m 


12.0 


DFLl 


0 


lndole-3-acetic acid amido synthetase 


6a 


ppa003113m 


3.8 


EIL3 


0 


Ethylene-insensitive3-like 3 




ppa001557m 


6.8 


ARF4 


0 


Auxin response factor 4 




ppa002082m 


^,1 


ARFIO 


0 


Auxin response factor 10 


6b 


ppa009583m 


8.7 


HDT3 


3x10-** 


Histone deacetylase 2C 




ppa006053m 


9.1 


HDA9 


0 


Histone deacetylase 9 




ppa007108m 


11.6 


HDA8 


0 


Histone deacetylase 8 




ppa022266m 


13.7 


AREB3 


2x10-=' 


ABA-responsive element binding protein 3 


7b 


ppa024294m 


12.4 


ATXR4 


4x10-"° 


Histone-lysine N-methyltransferase 




ppa026273m 


13.4 


At5g42560 


1x10-'" 


ABA-responsive protein {TB2/DP1, HVA22) 




ppaOl 2369m 


14.3 


TFLl 


7x10-"^ 


Phosphatidylethanolamine-binding protein 




ppaOl 5643m 


14.5 


CBF2 


1xlO-"8 


C-repeat/DRE binding factor 2 




ppa009356m 


14.7 


HDT3 


6x10-=' 


Histone deacetylase 2C 




ppa001213m 


15.4 


CLE 


0 


Histone-lysine N-methyltransferase 




ppa002248m 


15.6 


ABAl 


0 


Zeaxanthin epoxidase 




ppaOl 2239m 


16.4 


At4g37220 


7x10-^= 


Cold acclimation protein WCOR413 family 




ppa001803m 


16.7 


DDMl 


0 


Chromatin remodelling factor 
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Table 6 Putative orthologs in peach of Arabidopsis genes involved in chromatin modification and flowering located in 
QTL regions, identified by reciprocal blast analysis 

Complex or pathway Arabidopsis protein Gene model Putative peach ortholog ^'^^*P Genomic location QTL region 



Polycomb (PRC2) 


CLF 


At2g23380 


ppa001213m 


0 


scaffold_^ 


7: 15427185 


- 1 5434449 


7b 




VI L2 


At4g30200 


ppa001943m 


- ^-1 77 

6x10 


scaffold_' 


1: 5841849 - 


5849018 


la 


Polycomb (PRC1) 


EMF1 


At5g 11530 


ppa000318m 


2x10"^^ 


scaffold_' 


1: 37002560 


- 37008224 


lb 


Trithorax 


PKL 


At2g25170 


ppa000228m 


0 


scaffold_' 


1: 33913537 


- 33925777 


lb 




ASH2R 


At1g51450 


ppa006686m 


2x10"^^^ 


scaffold_' 


1: 1 940241 3 


- 1 9406234 


4c 


JMJ domain 


REF6 


At3g48430 


ppa000214m 


0 


scaffoldj 


y. 14087553 


- 14096069 


5b 


Histone acetylation 


HAC1 


At1g79000 


ppa000162m 


0 


scaffold_; 


?: 1240769 - 


1 249807 


3a 




HAM2 


At5g09740 


ppa005747m 


0 


scaffold_' 


1: 40838319 


- 40843052 


lb 




HDA6 


At5g63110 


ppa005308m 


0 


scaffoldj 


y. 14461481 


- 14464755 


5b 


Flowering 


SVP 


At2g22540 


ppaOl 1063m 


8x10^*^ 


scaffoldj 


5: 18701564 


- 1 8704722 


6c 




LFY 


At5g61850 


ppa006372m 


9x10"**' 


scaffold.! 


y. 13116676 


- 13119151 


5b 




TFL1 


At5g03840 


ppaOl 2369m 


1x10"'^ 


scaffoldj 


7; 14280983 


- 14282266 


7b 



dormancy release (EnD and EcD), and blooming date 
(BD) variables constitute a first group with high correl- 
ation values and similar QTLs. Additionally, the fact that a 
significant QTL for fruit harvest date also co-localizes with 
the major QTL for chilling requirement in LGl indicates 
that chilling requirement is a primary determinant of the 
reproductive phenology in peach. Heat requirement HREc 
and HRB, and PEnEc and PEnB form a second group of 
traits with common features. Finally, HD trait deserves a 
separate discussion due to the particular contribution of 
fruit developmental programs to this trait, in spite of its 
significant correlation with chilling requirement trait. 

Candidate genes for chilling requirement, dormancy 
release and blooming time 

QTLs for this first group of traits clustered in seven differ- 
ent map zones that corresponded to seven genomic regions 
defined by SNPs (lb, 3a, 3b, 4a, 6b, 6c and 7b). Some of 
these regions were identified in previous QTL works in Pm- 
nus species. Preceding attempts to describe loci affecting 
the blooming date trait in peach found QTLs co-localizing 
with our QTLs in regions lb [16] and 7b [13,19], whereas 
more recent studies dissecting the chilling requirement and 
blooming time traits in apricot, peach and almond identi- 
fied QTLs overlapping with our QTLs in regions lb, 4a, 6b 
and 7b [22-24]. The QTL in LG7 appeared consistently in 
most of these reports. Our analysis pointed to the presence 
of two adjacent QTL clusters that could contribute to both 
chilling requirement and dormancy-blooming time traits in 
this zone (Table 4). Due to their shortly overlapping posi- 
tions, these two putative clusters were fused and considered 
as one single QTL in the subsequent analysis (Figure 3). In 
almond, two QTLs were also found in adjacent positions in 
LG7 [24]. By assuming a high degree of synteny between al- 
mond and peach, the first QTL for flowering time identified 



in LG7 in almond co-localized with our QTL in region 7b, 
whereas the second QTL identified by these authors associ- 
ated to chilling and heat requirement traits did not overlap 
with region 7b. 

The availability of the peach genome sequence [34] fa- 
cilitates the identification of candidate genes by in silico 
search of genes within QTL intervals. Following this ap- 
proach, genes involved in light signalling, circadian 
clock, flowering regulation, cell cycle and phytohormone 
response were previously identified as candidate genes 
for bud phenology traits in poplar and apple [9-12]. 

Recent studies in apricot, peach and almond proposed 
DAM genes within the evergrowing [evg) locus as the 
most promising candidate genes for the major QTL af- 
fecting chilling requirement and blooming time in LGl 
[22-24], based on the genomic location of evg and the 
abundant literature conferring DAM genes a relevant 
role in bud dormancy maintenance. DAMl-6 genes are a 
set of six tandemly repeated MADS-box genes related to 
SHORT VEGETATIVE PHASE {SVP) of Arabidopsis 
thaliana that have been found partially deleted in the 
evg peach mutant showing non-dormant behaviour 
[25,35]. DAM genes are specifically expressed in buds 
and are affected differently by photoperiod and chilling 
signals [36]. DAMS and DAM6 expression correlated 
with the dormancy state of buds, being higher in dor- 
mant buds and lower after the fulfilment of chilling re- 
quirements prior to dormancy release [37-39]. The 
expression of DAMl, DAMS and DAM6 is also re- 
pressed during chilling stratification of the embryo, sug- 
gesting their participation in seed dormancy release 
mechanisms [40]. At the functional level, the heterol- 
ogous expression of DAMl gene from leafy spurge {Eu- 
phorbia esula) delayed flowering in Arabidopsis [41], 
and PmDAM6 from Japanese apricot {Prunus mume) led 
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to growth cessation and bud set in poplar under envir- 
onmental conditions favourable for growth [42]. 

Other candidates to be the major determinant of bud 
phenology located in genomic region lb are listed in 
Table 5. Among them, we have proposed chromatin re- 
modelling and modification factors such as PICKLE-like 
(ppa000228m), a putative SWI3C-like element of the 
SWITCH/SUCROSE NONFERMENTING (SWI/SNF) 
remodelling complex (ppa001566m), an HDA2-like his- 
tone deacetylase (ppa006590m), HAM2-like histone 
acetyltransferase (ppa005747m), an EARLY FLOWER- 
ING IN SHORT DAYS (EFS)-like histone methyltrans- 
ferase (ppa000056m), and EMBRYONIC FLOWERl 
(EMFl)-like component of the Polycomb Repressive 
Complexl (PRCl) (ppa000318m). The vernalization 
mechanisms converging on the expression of FLOWER- 
ING LOCUS C (FLC) in Arabidopsis have been pro- 
posed to share molecular features with the chilling- 
dependent release of bud dormancy mediated by DAM 
genes [43,44]. The chromatin modification mechanisms 
involved in FLC regulation include synthesis of non- 
coding RNAs, histone acetylation, trimethylation of 
H3K4, methylation of H3K36 by EES, trimethylation of 
H3K27 by PRC2 complex, and monoubiquitination of 
H2A by PRCl among others [32,45]. Interestingly 
DAMl from leafy spurge and DAM6 from peach are 
regulated at the chromatin level by demethylation of 
H3K4 and trimethylation of H3K27 following chilling 
accumulation, in a similar way to FLC [30,41]. In 
addition DAM6 chromatin also showed chilling- 
dependent differences in H3 acetylation [30]. Altogether, 
these and other published data in chestnut [46,47] 
emphasize a prominent role of chromatin modifying path- 
ways in bud dormancy mechanisms. 

Other candidate genes in region lb are putative 
components of the ABA signalling pathway, such as 
ppa003748m coding for a protein phosphatase 2C (PP2C) 
simUar to HYPERSENSITIVE TO ABAl (HABl), and 
ppa006503m and ppa024363m coding for proteins simi- 
lar to ABA-RESPONSIVE ELEMENT BINDING PRO- 
TEIN 1 (AREBl) and AREB3 [48]. Additional HABl- 
like and A7?£B3-like genes are also found in regions 3a 
and 6b. HABl and other related PP2Cs perform a cen- 
tral role in the negative regulation of ABA signalling in 
Arabidopsis, which is overcome by the ABA-dependent 
interaction of PP2Cs with the ABA-receptor PYL5 [49]. 
In contrast to the well-established role of ABA in seed 
dormancy processes, only few molecular data support 
the function of ABA in promoting and maintaining dor- 
mancy in buds [2,3,50,51]. Furthermore, manipulating 
the expression of the poplar ortholog of ABSCISIC 
ACID INSENSITIVE 3 {ABI3) caused alterations in bud 
formation and misregulation of numerous genes in buds 
[52]. 



In region lb we have also identified putative 
flowering-related genes, such as ppa013757m similar to 
FLOWERING PROMOTING FACTOR 1 (FPFl) [53], 
and putative regulatory genes found up-regulated in la- 
tent buds such as ppa008979m and ppa012329m [30]. 

Genomic regions 3a, 6b and 7b were also considered 
important for the chilling requirement trait. Among the 
candidate genes present in these regions, we found other 
chromatin-related factors, such as ppa001213m, the 
peach ortholog of CURLY LEAF (CLF), a component of 
the PRC2 complex involved in the trimethylation of his- 
tone H3 at lysine 27 [54]. As already proposed in previ- 
ous works, PRC2 complexes could contribute to bud 
dormancy release in Prunus species [55,56], and more 
specifically to H3K27 trimethylation observed in DAM6 
concomitantly with gene down-regulation [30]. The 
genes ppa004975m and ppa024294m codify for other 
putative histone methyltransferases containing the SET- 
domain, with similarity to Arabidopsis SDG40 and 
ATXR4 respectively. In region 7b, ppa001803m codes for 
a putative SWI2/SNF2 chromatin-remodelling ATPase 
similar to DDMl, which makes the heterochromatin 
bound to histone HI accessible to DNA methyltransfer- 
ases [57]. 

Further histone acetyltransferases and deacetylases are 
localized in regions 3a, 6b and 7b. Among them, 
ppa009583m and ppa009356m show similarity to HDT3 
gene, coding for a histone deacetylase that modulates 
the ABA response [58]. Other putative elements of 
ABA signalling, ABA biosynthesis and stress response 
pathways in peach are ppa004252m, ppa002515m, 
ppa022266m, ppa026273m, ppa015643m, ppa002248m 
and ppaO 12239m. 

Candidate genes for heat requirement 

In contrast to the well-established genetic component of 
chilling requirements and the close relationship between 
flowering time and chilling requirements described so 
far, the genetic control of heat requirements in Prunus 
species has been a matter of discussion in the bud dor- 
mancy field. Couvillon & Erez [59] considered the varia- 
tions in heat requirement to be due to excessive chilling 
and found no genetic differences in heat requirements 
among cultivars. However, the negative correlation 
found previously between chilling and heat requirements 
has been argued to suggest the existence of a potential 
contribution of genetic factors to the heat requirement 
trait [23,60]. We have observed a similar negative correl- 
ation of PEnEc and PEnB with chilling requirement vari- 
ables in this work; however HREc and HRB traits were 
not related significantly to chilling requirements with 
the exception of a positive correlation found between 
HREc and CRW/CRU/CRD in AA location (Table 2). 
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We found 13 significant QTLs for HREc, HRB, PEnEc 
and PEnB, located in nine different genomic regions 
(Table 4, Figure 3). Seven QTLs overlapped with chilling 
requirement QTLs in regions lb, 3a, 6b and 7b. No co- 
incidences with previous reports were observed, with the 
exception of QTLs in the genomic region lb [23]. 

In addition to candidate genes proposed for genomic 
regions containing chilling requirement QTLs, commen- 
ted in the previous section, a reciprocal blast analysis for 
the search of peach genes orthologous to chromatin and 
flowering genes from Ambidopsis resulted in the candi- 
date gene list presented in Table 6. The transcript model 
ppa001943m, located in region la, was a putative ortho- 
log of VERNALIZATION INSENSITIVE 3-LIKE 2 (VIL2), 
coding for a component of PRC2 complexes involved in 
flowering under non-inductive conditions through the 
epigenetic regulation of the floral repressor MADS AF- 
FECTING FLOWERING 5 {MAF5) [61]. In region 4c we 
identified a putative ortholog of ARABIDOPSIS ASH2 
RELATIVE (ASH2R) gene, a regulator of flowering time 
required for H3K4 trimethylation and the proper expres- 
sion of FLC and FLC homologs [62]. Finally, in region 
5b we identified the putative orthologs of the chromatin 
regulators RELATIVE OF EARLY FLOWERING 6 (REF6) 
and HDA6 [63,64], and the floral modulator LEAFY 
(LEY) [65]. 

Candidate genes for fruit harvest date 

The analysis of the harvest date trait resulted in a major 
QTL in region 4b showing numerous precedents in re- 
lated works in peach and apricot [16,19,20]. The tran- 
script model ppa010982m, similar to ETHYLENE 
RESPONSIVE ELEMENT BINDING FACTOR 4 {ERF4) 
from Ambidopsis, has been already proposed as a candi- 
date gene for fruit maturation date in different Prunus 
species [20]. ERF4 is a transcriptional repressor modulat- 
ing ethylene and ABA responses in Arabidopsis [66]. In 
peach different ERF genes have been found up-regulated 
in ripening fruit [67], whereas similar ERFs have been 
postulated to be involved in fruit ripening regulation in 
apple [68]. 

On the same region 4b, the gene ppa022739m codes 
for a putative transcription factor containing the 
Squamosa-Promoter Binding Protein (SBP) domain, 
present in the tomato fruit ripening factor COLORLESS 
NON-RIPENING (CNR) [69]. 

However, recent fine mapping of a locus controlling 
maturity date in two segregating populations of peach 
limited the search to a 220 kb stretch within region 4b. 
The maturity date locus co-segregated with an indel into 
the gene ppa008301m coding for a NAC type transcrip- 
tion factor, which points to this gene as a firm candidate 
for controlling ripening time in peach [70]. 



In region 6a we should emphasize the presence of 
ppa003113m gene, with similarity to ETHYLENE-IN- 
SENSITIVE3-LIKE 3 {EIL3), involved in regulation of the 
sulfur-limitation response in Ambidopsis [71] and similar 
to elements of the ethylene pathway. Other candidate 
genes in regions 4b and 6a were hypothetically related to 
auxin synthesis and response (ppa002986m, ppa001557m 
and ppa002082m), since auxin is known to be involved in 
fruit set and ripening [72]. 

Conclusions 

This work was aimed at the identification of genetic fac- 
tors conditioning the phenological behaviour of peach. 
We have identified QTLs for nine traits related to bud 
dormancy, flowering and fruit harvest in a hybrid popu- 
lation of peach in two different locations. QTLs were lo- 
cated in a SNP-based genetic linkage map. A search of 
candidate genes for these QTLs rendered different genes 
related to flowering regulation, chromatin modification 
and hormone signalling. Additional studies including the 
characterization of proposed candidate genes in germ- 
plasm collections and functional approaches are re- 
quired to identify the genes involved in dormancy, 
blooming and fruit maturation among these lists of can- 
didate genes. The characterization of natural alleles of 
these genes might offer molecular tools to predict the 
potential performance of different Prunus species and 
cultivars under changing climatic conditions. 

Methods 

Plant material 

The plant material used in the study was a progeny of 
107 individuals derived from a cross carried out during 
2008 between the Fl selection 'V6' (named MxR Ol in 
[73]) and the Brazilian non-melting peach cultivar 
'Granada'. 'V6' selection was derived from a cross per- 
formed in 2005 between the Spanish non-melting peach 
cultivar 'Maruja' (high chilling requirement and late rip- 
ening) and the North- American melting peach cultivar 
'Red Candem' (low chilling and early ripening). Two in- 
dividuals per genotype were grafted on 'Garnem' (hybrid 
almond x peach) rootstock, and then cultivated in the 
experimental orchards Agua Amarga (AA; 38° 18' 41" N 
1° 31' 31" W, 344 m over sea level) and El Jimenado (EJ; 
37° 45' 31" N 1° 01' 35" W, 80 m over sea level), both of 
them situated in the Region of Murcia, at the southeast 
of Spain. The experimental orchards represented differ- 
ent climatic conditions regarding chilling and heat 
accumulation, maximum, minimum and medium tem- 
peratures and humidity during winter and spring. From 
the 107 progeny trees 86 genotypes were grown at EJ, 74 
at AA, and 70 genotypes were common to both loca- 
tions. Horticultural practices such as pruning, irrigation, 



Romeu ef al. BMC Plant Biology 2014, 14:52 
http://www.biomedcentral.com/1471-2229/14/52 



Page 13 of 16 



fertilization and control of weeds, insects and diseases 
were consistently performed at both orchards. 

Phenotypic assessment 

Chilling requirements, heat requirements and blooming 
dates were evaluated in the winter and spring of 2012; 
fruit harvest dates were measured in the spring of 2011 
and 2012. To assess chilling requirements for endodor- 
mancy release, nine one-year old shoots per genotype 
with a length of 20-25 cm were picked weekly. Groups 
of three shoots were placed in bottles containing dis- 
tilled water with 3% sucrose, and incubated in a growth 
chamber subject to 12 h photoperiod at 22°C. The basal 
end of shoots was cut and the water renewed once per 
week. We considered that endodormancy was completed 
when the percentage of buds that reached the green 
stage (stage C) according to the Baggiolini code [74] was 
higher than 30% in the three groups of shoots after 
10 days. It was difficult to obtain higher percentages of 
bud break under the artificial conditions (cutting shoots) 
of the laboratory. Quantification of chilling accumula- 
tion at the dormancy release date was performed using 
the three most common models: Weinberger [75], Utah 
[76] and dynamic model [77,78]. Hourly air tempera- 
tures were recorded from beginning of winter to harvest 
date by the SIAM station [79] and by temperature sen- 
sors (Testo T174). 

To evaluate ecodormancy release and blooming time, 
the main phenological stages [74] were visually identified 
on the field weekly, from the beginning of winter to fruit 
set in 2011-2012. The different phenological stages were 
assessed quantitatively based on the ratio of buds. Thus, 
we considered that ecodormancy was released when 50% 
of buds had reached the green stage (stage C) in the 
field, and a percentage of 50% of open flowers (stage F) 
served to establish blooming date. Both dates were fur- 
ther expressed as Julian days and periods between endo- 
dormancy, ecodormancy and blooming time were 
calculated. 

Heat requirements were calculated as the growing de- 
gree hours (GDH) accumulated from the release of 
endodormancy to the ecodormancy release and bloom- 
ing dates following the Anderson model [80]. The har- 
vest date was determined in situ based on fruit colour 
and firmness. 

Statistical analyses 

Statistical analyses were performed using the Statgraphics 
5.1 package (Statpoint Technologies, Warrenton, VA, 
USA). All correlations between traits were calculated 
using the Pearson coefficient. Correlations between traits 
in the EJ and AA locations employed exclusively those ge- 
notypes present in both locations. Departure from the 



normal distribution of traits was assessed by the calcula- 
tion of skewness and kurtosis of frequency distributions. 

SNP genotyping and map construction 

DNA from the parentals and progeny were extracted 
from 50 mg of leaf tissue following the method of Doyle 
& Doyle [81]. The concentration of DNA was checked 
by comparison with standard DNA ladders in agarose 
gels and with Quant-iTTM PicoGreen H Assay (Life 
Technologies, Grand Island, NY, USA). Samples were 
genotyped using the International Peach SNP Consor- 
tium (IPSC) peach 9 K Infinium" II array [27] at the 
Genotyping and Genetic Diagnosis Unit (Health Re- 
search Institute, INCLIVA, Valencia, Spain). The SNP 
array includes information about the physical position of 
all the SNPs in 9 genome scaffolds, being the first eight 
ones corresponding to the eight peach chromosomes. 
After visual inspection of genotype calls, monomorphic 
SNPs and SNPs with more than 5% of missing data were 
removed. The map construction has been described in 
[28] and it will be published with further details else- 
where. Briefly, we followed the two-way pseudo-test 
cross approach [82]. Homozygous SNPs in one parent 
and heterozygous in the other parent were selected to 
generate a genetic map for each parent, discarding SNPs 
heterozygous for both parents, as these markers were 
not used for QTL mapping and an integrated map was 
not necessary because the physical position of the SNPs 
was already known. A total of 1,970 SNPs segregated 
(1:1) for the 'V6' parent and 895 for 'Granada'. From this 
data set, we removed the SNPs that showed exactly the 
same genotypic segregation to obtain a non-redundant 
and simplified map more suitable for QTL mapping. 
Marker data was coded as cross-pollinator (CP) and 
linkage analysis was performed with JoinMap" 4 [83] 
with a minimum LOD from 6.0 to 8.0. Map construction 
was performed using the regression mapping algorithm 
[83] and the default JoinMap' parameters (Rec = 0.40, 
LOD = 1, Jump = 5.0, and ripple = 1). The genotyping data 
was re-coded as pseudo back-cross and the order of the 
markers was double checked with Mapmaker 3.0 [84]. 
The Kosambi mapping function was used to convert re- 
combination frequencies into map distances. The maps 
for each parent were drawn with MapChart 2.2 [85]. 

QTL analysis 

In order to facilitate computer analysis, the genetic linkage 
was condensed, eliminating SNPs that mapped in the 
same position or very close (i.e. less than 2 cM). The maps 
from each parent were analyzed independently and coded 
as two independent backcross populations. QTL analysis 
was performed with WinQTLcartographer 2.5 [86] by 
Composite Interval Mapping (CIM), and the LOD thresh- 
old to declare a QTL significant at P < 0.05 was calculated 
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by a permutation test for the whole genome [87], imple- 
mented in WinQTLcartographer. A two-LOD support 
interval was taken as a confidence interval for the detected 
QTLs. That is, the confidence intervals were limited by a 
decrease of two in the LOD score at both sides of the 
QTL peal<. 

Candidate gene selection 

By using the BioMart tool in the phytozome web-page 
[88] we obtained the annotated transcript models con- 
tained between the SNP markers flanking the QTL re- 
gions (Additional file 1 and Additional file 2). These 
markers included the two-LOD confidence interval of 
their respective QTL. Those genes similar to known 
genes involved in ABA regulation, cold acclimation, 
ethylene signalling, chromatin modification, flowering 
and vernalization pathways within the most relevant re- 
gions lb, 3a, 4b, 6a, 6b and 7b were selected as candi- 
date genes. 

In order to identify putative orthologs in peach of Ara- 
bidopsis genes related to chromatin modification and 
flowering pathways we performed a reciprocal blast ana- 
lysis at phytozome [88]. First we made a blastp similarity 
search by using the protein sequence of selected genes 
as query. The first hit in the peach genome was subse- 
quently compared with the Arabidopsis genome by 
blastp search, and those genes found reciprocally by the 
searches in both the peach and Arabidopsis genomes 
were considered as putative orthologs. 

Additional files 



Additional file 1: Name and genomic position of SNPs flanking the 
QTL regions. 

Additional file 2: Annotation of genes within the QTL regions. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

JFR carried out the genetic studies and helped to draft the manuscript. AJM 
carried out the linkage map construction and QTL analysis. GS participated in 
the linkage map construction and QTL analysis. AG participated in the 
design and coordination of the study. JGB participated in the design and 
coordination of the study. MLB conceived the study, and participated in its 
design and coordination. GR carried out the candidate gene search and 
drafted the manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

We thank Matilde Gonzalez for technical assistance. This work was supported 
by the Instituto Nacional de Investigacion y Tecnologia Agraria y Alimentaria 
(INIA)-FEDER (grant no. RTA2007-00060), and the Ministry of Science and 
Innovation of Spain (grant no. AGL201 0-20595). 

Author details 

'instituto Murciano de Investigacion y Desarrollo Agrario y Alimentario 
(IMIDA), La Alberca, Murcia, Spain. ^Institute de Biologia Molecular y Celular 
de Plantas (IBMCP, CSIC-UPV), Valencia, Spain. ^Instituto Valenciano de 
Investigaciones Agrarias (MA), Moncada, Valencia, Spain. 



Received: 10 December 2013 Accepted: 18 February 2014 
Published: 22 February 2014 

References 

1. Lang GA: Dormancy: a new universal terminology. Hon Sci 1987, 
22:817-820. 

2. Rohde A, Bhalerao RP: Plant dormancy in the perennial context. 
Trends Plant Sci 2007, 12:217-223. 

3. Cooke IFK, Friksson ME, Junttila 0: The dynamic nature of bud dormancy 
in trees: environmental control and molecular mechanisms. Plant Cell Env 
2012, 35:1707-1728. 

4. Coville FV: The influence of cold in stimulating the growth of plants. 
Proc Natl Acad Sci USA] 920, 6:434-435. 

5. Chuine I: Why does phenology drive species distribution? Phllos Trans 
R Soc Land B Biol Sci 2010 365:3149-3160 

6. Hanninen H, Tanino K: Tree seasonality in a warming climate. Trends Plant 
Sd 2011, 16:412-416. 

7. Yu H, Luedeling E, Xu J: Winter and spring warming result in delayed 
spring phenology on the Tibetan Plateau. Proc Natl Acad Sci U S A20]0, 
107:22151-22156 

8. Nicotra AB, Atkin OK, Bonser SP, Davidson AlVl, Finnegan EJ, Mathesius U, 
Poot P Purugganan MD, Richards CU Valladares F, van Kleunen M: Plant 
phenotypic plasticity in a changing climate. Trends Plant Sci 2010, 
15:684-692. 

9 Frewen BE Chen THH, Howe GT Davis J, Rohde A, Boerjan W, Bradshaw HD 
Jr: Quantitative trait loci and candidate gene mapping of bud set and 
bud flush in Populus. Genetics 2000, 154:837-845. 

1 0. Rohde A, Storme V, lorge V, Gaudet M, Vitacolonna N, Fabbrini F, Ruttink T 
Zaina G, Marron N, Dillen S Steenackers M, Sabatti M, Morgante M, Boerjan 
W, Bastien C: Bud set in poplar - genetic dissection of a complex trait in 
natural and hybrid populations. New Phytol 201 1, 189:106-121. 

1 1 . Fabbrini F, Gaudet M, Bastien C, Zaina G, Harfouche A, Beritognolo I, Marron 
N, Morgante M, Scarascia-Mugnozza G, Sabatti M: Phenotypic plasticity, 
QTL mapping and genomic characterization of bud set in black poplar. 
BMC Plant Biol 2012, 12:47. 

12. Celton J-M, Martinez S, Jammes M-J, Bechti A, Salvi S, Legave 1-M, Costes E: 
Deciphering the genetic determinism of bud phenology in apple 
progenies: a new insight into chilling and heat requirement effects on 
flowering dates and positional candidate genes. New Phytol 201 1, 
192:378-392. 

13. Diriewanger E, Moing A, Rothan C, Svanella L, Pronier V, Guye A, Plomion C, 
Monet R: Mapping QTL controlling fruit quality in peach (Prur)us persica 
(L) Batsch). TheorAppI Genet 1999 98:18-31. 

14. Ballester J, Socias i Company R, Arus P, de Vicente MC: Genetic mapping of 
a major gene delaying blooming time in almond. Plant Breed 2001, 
120:268-270 

15. Verde I, Quarta R Cedrola C, Dettori MT: QTL analysis of agronomic traits 
in a BC, peach population. Acta Hort 2002, 592:291-297. 

16. Quilot B, Wu BH, Kervella 1, Genard M, Foulongne M, Moreau K: QTL 
analysis of quality traits in an advanced backcross between Prur]us 
persica cultivars and the wild relative species P. davidiana. Theor AppI 
Genet 2004 109:884-897. 

17. Sanchez-Perez R, Howad W, Dicenta F, Arus P, Martinez-Gomez P: Mapping 
major genes and quantitative trait loci controlling agronomic traits in 
almond. Plant Breed 2007, 126:310-318. 

18. Campoy JA, Ruiz D, Egea J, Rees DIG, Celton JM, Martinez-Gomez P: 
Inheritance of flowering time in apricot {Prunus armeniaca L.) and 
analysis of linked quantitative trait loci (QTL) using simple sequence 
repeat (SSR) markers. Plant IVIol Biol Rep 2011, 29:404-410 

19. Eduardo I, Pacheco I, Chietera G, Bassi D, Pozzi C, Vecchietti A, Rossini L 
QTL analysis of fruit quality traits in two peach intraspecific populations 
and importance of maturity date pleiotropic effect. Tree Genet Genomes 
2011, 7:323-335 

20. Diriewanger E, Quero-Garcia J, Le Dantec L, Lambert P, Ruiz D, Dondini L, 
lla E, Quilot-Turion B, Audergon JM, Tartarini S, Letourmy P, Arus P: 
Comparison of the genetic determinism of two key phenological traits, 
flowering and maturity dates, in three Prunus species: peach, apricot and 
sweet cherry. Heredity imi, 109:280-292. 

21. Socquet-Juglard D, Christen D, Devenes G, Gessler C, Duff/ B, Patocchi A: 
Mapping architectural, phenological, and fruit quality QTLs in apricot. 
Plant IVIol Biol Rep 2013, 31:387-397. 



Romeu ef al. BMC Plant Biology 2014, 14:52 
http://www.biomedcentral.com/1471-2229/14/52 



22. Olukolu BA, Trainin T, Fan S, Kole C, Bielenberg DG, Reighard GL, Abbott AG, 
Holland D: Genetic linl<age mapping for molecular dissection of chilling 
requirement and budbreak in apricot [Prunus armeniaca L). Genome 

2009, 52:819-828. 

23. Fan S, Bielenberg DG, Zhebentyayeva TN, Reighard GL, Okie WR, Holland D, 
Abbott AG: IVlapping quantitative trait loci associated with chilling 
requirement, heat requirement and bloom date in peach {Prunus 
persica). New Phytol 2010, 185:917-930. 

24. Sanchez-Perez R, Dicenta F, Martinez-Gomez P: Inheritance of chilling and 
heat requirements for flowering in almond and QTL analysis. Tree Genet 
Genomes 2012, 8:379-389. 

25. Bielenberg DG, Wang Y, Li Z, Zhebentyayeva T, Fan S, Reighard GL, Scorza R, 
Abbott AG: Sequencing and annotation of the evergrowmg locus in 
peach (Prur)us persica [L] Batsch) reveals a cluster of six MADS-box 
transcription factors as candidate genes for regulation of terminal bud 
formation. Tree Genet Genomes 2008, 4:495-507. 

26. Jimenez S, Li Z, Reighard GL, Bielenberg DG: Identification of genes 
associated with growth cessation and bud dormancy entrance using a 
dormancy-incapable tree mutant. BMC Plant Biol 2010, 10:25. 

27. Verde I, Bassil N, Scalabrin S, Gilmore B, Lawley CT, Gasic K, Micheletti D, 
Rosyara UR, Cattonaro F, Vendramin E, IVlain D, Aramini V, Bias AL, Mocker 
TC, Bryan DW, Wilhelm L, Troggio M, Sosinski B, Aranzana MJ, Arus P, lezzoni 
A, Morgante M, Peace C: Development and evaluation of a 9 K SNP array 
for peach by internationally coordinated SNP detection and validation in 
breeding germplasm. PLoS One 2012, 7:e35668. 

28. Sanchez G: Cartografiado de QTLs y Genes Candidatos Asociados a Metabolitos 
Determinantes de la Calidad de Fmto en Melocotdn, PhD Thesis. Universidad 
Politecnicade Valencia: 201 3. http://hdl.handle.net/10251/34511. 

29. Leida C, Terol J, Marti G, Agustf M, Llacer G, Badenes ML, Rios G: Identification 
of genes associated with bud dormancy release in Prunus persica by 
suppression subtractive hybridization. Tree Physiol 20]0, 30:655-566. 

30. Leida C, Conesa A, Llacer G, Badenes ML, Rios G: Histone modifications 
and expression of DAM6 gene in peach are modulated during bud 
dormancy release in a cultivar-dependent manner. New Phytol 2012, 
193:67-80. 

31. Leida C, Romeu JF, Garcfa-Brunton J, Rios G, Badenes ML: Gene expression 
analysis of chilling requirements for flower bud break in peach. Plant 
Breed 201 2, 131:329-334. 

32. Holec S, Berger F: Polycomb group complexes mediate developmental 
transitions in plants. Plant Physiol 201 2, 1 58:35-43. 

33. Pandey R, Muller A, Napoli CA, Selinger DA, Pikaard CS, Richards FJ, Bender 
J, Mount DW, Jorgensen RA: Analysis of histone acetyltransferase and 
histone deacetylase families of Arabidopsis thaliana suggests functional 
diversification of chromatin modification among multicellular 
eukaryotes. Nucleic Adds Res 2002, 30:5036-5055. 

34. The International Peach Genome Initiative: The high-quality draft genome 
of peach {Prunus persica) identifies unique patterns of genetic diversity, 
domestication and genome evolution. Nature Genet 2013, 45:487-494. 

35. Jimenez S, Lawton-Rauh AL, Reighard GL, Abbott AG, Bielenberg DG: 
Phylogenetic analysis and molecular evolution of the dormancy 
associated MADS-box genes from peach. BMC Plant Biol 2009, 9:81. 

36 Li Z, Reighard GL, Abbott AG, Bielenberg DG: Dormancy-associated MADS 
genes from the EVG locus of peach {Prunus persica [L] Batsch) have 
distinct seasonal and photoperiodic expression patterns. J Exp Bot 2009, 
60:3521-3530. 

37. Yamane H, Ooka T, Jotatsu H, Sasaki R, Tao R: Expression analysis of 
PpDAMS and PpDAIV\6 during flower bud development in peach 
{Prunus persica). Sci Hortic 201 1, 129:844-848. 

38. Jimenez S, Reighard GL, Bielenberg DG: Gene expression of DAMS and 
D/IM6 is suppressed by chilling temperatures and inversely correlated 
with bud break rate. Plant Mol Biol 2010, 73:157-167. 

39. Yamane H, Ooka T, Jotatsu H, Hosaka Y, Sasaki R, Tao R: Expressional 
regulation of PpDAMS and PpDAM6, peach {Prunus persica) 
dormancy-associated MADS-box genes, by low temperature 

and dormancy-breaking reagent treatment. J Exp Bot 201 1, 62:3481-3488. 

40. Leida C, Conejero A, Arbona V, Gomez-Cadenas A, Llacer G, Badenes ML, Rios G: 
Chilling-dependent release of seed and bud dormancy in peach associates 
to common changes in gene expression. PLoS One 201 2, 7:e35777. 

41. Horvath DP, Sung S, Kim D, Chao W, Anderson J: Characterization, 
expression and function of DORMANCY ASSOCIATED MADS-BOX genes 
from leafy spurge. Plant Mol Biol 2010, 73:169-179. 



Page 15 of 16 



42. Sasaki R, Yamane H, Ooka T, Jotatsu H, Kitamura Y, Akagi T, Tao R: 
Functional and expressional analyses of PmDAM genes associated with 
endodormancy in Japanese apricot. Plant Physiol 201 1, 157:485-497. 

43. Horvath D: Common mechanisms regulate flowering and dormancy. 
P/onf Sd 2009, 177:523-531. 

44. Hemming MN, Trevaskis B: Make hay when the sun shines: the role of 
MADS-box genes in temperature-dependent seasonal flowering 
responses. Plant Sci 201 1 , 1 80:447-453. 

45. He Y: Chromatin regulation of flowering. Trends Plant Sd 201 2, 
17:556-562. 

46. Santamarfa MF, Hasbun R, Valera MJ, Meijon M, Valledor L, Rodriguez JL, 
Toorop PF, Canal MJ, Rodriguez R: Acetylated H4 histone and genomic 
DNA methylation patterns during bud set and bud burst in Castanea 
sativa. J Plant Physiol 2009, 1 66: 1 360- 1 369. 

47. Santamaria MF, Rodriguez R, Canal MJ, Toorop PE: Transcriptome analysis 
of chestnut {Castanea sativa) tree buds suggests a putative role for 
epigenetic control of bud dormancy. Ann Bot 201 1, 108:485-498. 

48. Cutler SR, Rodriguez PL, Finkelstein RR, Abrams SR: Abscisic acid: 
emergence of a core signaling network. Annu Rev Plant Biol 2010, 
61:651-679 

49. Santiago J, Rodrigues A, Saez A, Rubio S, Antoni R, Dupeux F, Park SY, 
Marquez JA, Cutler SR, Rodriguez PL: Modulation of drought resistance by 
the abscisic acid receptor PYL5 through inhibition of clade A PP2Cs. 
Plant J 2009, 60:575-588. 

50. Arora R, Rowland LJ, Tanino K: Induction and release of bud dormancy in 
woody perennials: a science comes of age. HortSd 2003, 38:91 1 -921 . 

51. Horvath DP, Anderson JV, Chao WS, Foley ME Knowing when to grow: 
signals regulating bud dormancy. Trends Plant Sd 2003, 8:534-540. 

52. RuttinkT, Arend M, Morreel K, Storme V, Rombauts S, Fromm J, Bhalerao RP, 
Boerjan W, Rohde A: A molecular timetable for apical bud formation and 
dormancy induction in poplar. Plant Cell 2007, 19:2370-2390. 

53. Melzer S, Kampmann G, Chandler J, Apel K: FPFl modulates the 
competence to flowering in Arabidopsis. Plant J 1999, 18:395-405. 

54. Schubert D, Primavesi L, Bishopp A, Roberts G, Doonan J, Jenuwein T, 
Goodrich J: Silencing by plant Polycomb-group genes requires dispersed 
trimethylation of histone H3 at lysine 27. EMBO J 2006, 25:4638-4649. 

55. Fan S: Mapping Quantitative Trait Loci Associated with Chilling Requirement 
and Bloom Date in Peach, PhD Thesis. Clemson University: 2010. 
http://tigerprints.clemson.edu/all_dissertations/530. 

56. Olukolu BA: The Genetics of Chilling Requirement in Apricot (Prunus Armeniaca L), 
PhD Thesis. Clemson University: 2010. http://tigerprint5.clemson.edu/ 
all_dissertations/537. 

57. Zemach A, Kim MY Hsieh PH, Coleman-Derr D, Fshed-Williams L, Thao K, 
Harmer SL, Zilberman D: The Arabidopsis nucleosome remodeler DDM1 
allows DNA methyltransferases to access HI -containing heterochromatin. 
Ce//2013, 153:193-205. 

58. Sridha S, Wu K: Identification of AtHD2C as a novel regulator of abscisic 
acid responses in Arabidopsis. Plant J 2006, 46:124-133. 

59. Couvillon GA, Erez A: Influence of prolonged exposure to chilling 
temperatures on bud break and heat requirement for bloom of several 
fruit species. J Amer Soc Hort Sd 1 985, 1 1 047-50. 

60. Ruiz D, Campoy JA, Fgea J: Chilling and heat requirements of apricot 
cultivars for flowering. Environ Exp Bot 2007, 61:254-263. 

61 . Kim DH, Sung S: The Plant Homeo Domain finger protein, VIN3-LIKE 2, is 
necessary for photoperiod-mediated epigenetic regulation of the floral 
repressor, MAPS. Proc Natl Acad Sd U S A 20]0, ^ 07:1 7029-1 7034. 

62. Jiang D, Kong NC, Gu X, Li Z, He Y: Arabidopsis COMPASS-like complexes 
mediate histone H3 lysine-4 trimethylation to control floral transition 
and plant development. PLoS Genet 201 1 , 7:el 001 330. 

63. Lu F, Cui X, Zhang S, Jenuwein T, Cao X: Arabidopsis REF5 is a histone H3 
lysine 27 demethylase. Nat Genet 201 1 , 43:71 5-71 9. 

64 Yu CW, Liu X, Luo M, Chen C, Lin X, Tian G, Lu Q, Cui Y, Wu K HISTONE 
DEACETYLASE6 interacts with FLOWERING LOCUS D and regulates 
flowering in Arabidopsis. Plant Physiol 2011, 156:173-184. 

65. Weigel D, Alvarez J, Smyth DR, Yanofsky ME, Meyerowitz FM: LEAFY 
controls floral meristem identity in Arabidopsis. Cell 1992, 69:843-859. 

66. Yang Z, Tian L, Latoszek-Green M, Brown D, Wu K: Arabidopsis ERE4 is a 
transcriptional repressor capable of modulating ethylene and abscisic 
acid responses. Plant Mol Biol 2005, 58:585-596. 

67. Bonghi C, Trainotti L, Botton A, Tadiello A, Rasori A, Ziliotto F, Zaffalon V, 
Casadoro G, Ramina A: A microarray approach to identify genes involved 



Romeu ef al. BMC Plant Biology 2014, 14:52 Page 16 of 16 

http://www.biomedcentral.com/1471-2229/14/52 



69. 



70. 



72. 
73. 

74. 

75. 
76. 



79. 
80. 

81. 
82. 

83. 
84. 

85. 

86. 
87. 



in seed-pericarp cross-tall< and development in peach. BMC Plant Biol 
2011, 11:107. 

Wang A, Tan D, Takahashi A, Li 17, Harada T: IVldERFs, two ethylene-response 
factors involved in apple fruit ripening. J Exp Bot 2007, 58:3743-3748. 
Manning K, Tor M, Poole M, Hong Y, Thompson AJ, King GJ, Giovannoni JJ, 
Seymour GB: A naturally occurring epigenetic mutation in a gene 
encoding an SBP-box transcription factor inhibits tomato fruit ripening. 
Nat Genet 2006, 38:948-952. 

Pirona R, Eduardo i, Pacheco I, Da Silva Linge C, Miculan M, Verde I, Taitarini 
S, Dondini L, Pea G, Bassi D, Rossini L: Fine mapping and identification of 
a candidate gene for a major locus controlling maturity date in peach. 

BMC Plant Biol 2013, 13:166. 

Maruyama-Nakashita A, Nakamura Y, Tohge T, Saito K, Takahashi H: 
Arabidopsis SLIM1 is a central transcriptional regulator of plant sulfur 
response and metabolism. Plant Ceil 2006, 1 8:3235-325 1 . 
Seymour G, Poole M, Manning K, King GJ: Genetics and epigenetics of 
fruit development and ripening. Curr Opin Plant Biol 2008, 1 1 :58-63. 
Sanchez G, Besada C, Badenes ML, Monforte AJ, Graneli A: A non-targeted 
approach unravels the volatile network in peach fruit. PLoS One 2012, 
7:e38992. 

Baggiolini M: Stades reperes de la Abricotier- Stades Reperes de la 
Pecher. Stades Reperes du Ceresier - Stades reperes du Prunier. in Guide 
Pratique de Defense des Cultures. Paris: ACTA; 1 980. 
Weinberger JH: Chilling requirements of peach varieties. Proc Amer Soc 
HortSci 1950, 56:122-128. 

Richardson EA, Seeley SD, Walker DR: A model for estimating the 
completition of rest for "Redheaven" and "Elberta" peach. HortScI 1974, 
9:331-332. 

Fishman S, Erez A, Couviiion GA: The temperature dependence of 
dormancy breaking in plants - mathematical-analysis of a 2- step model 
involving a cooperative transition. J Theor Biol 1 987, 1 24:473-483. 
Fishman S, Erez A, Couviiion GA: The temperature dependence of 
dormancy breaking in plants - computer-simulation of processes studied 
under controlled temperatures. J Theor Biol 1987, 126:309-321. 
SIAIVl station. httpy/siam.imida.es/apex/f?p=l 01:1 10:657569513097788. 
Anderson JL, Richardson EA, Kesner CD: Validation of chill unit and flower 
bud phenology models for "Montmorency" sour cherry. Acta Hort 1986, 
184:71-78. 

Doyie JJ, Doyle JL: A rapid isolation procedure for small quantities of 
fresh leaf tissue. Phytochem Bull 1987, 19:1 1-15. 
Grattapagiia D, Sederoff R: Genetic-linkage maps of Eucalyptus grandis 
and Eucalyptus urophylla using a pseudo-testcrossmapping strategy and 
RAPD markers. Genetics 1994, 137:1 121 -11 37. 

Stam P: Construction of integrated genetic-linkage maps by means of a 
new computer package — Joinmap. Plant J 1993, 3:739-744. 
Lander ES, Green P, Abrahamson J, Barlow A, Daley M, Lincoln S, Newburg 
L: MAPIVIAKER: an interactive computer package for constructing primary 
genetic linkage maps of experimental and natural populations. Genomics 
1987, 1:174-181. 

Voorrips RE MapChart: software for the graphical presentation of linkage 
maps and QTLs. J Hered 2002, 93:77-78. 

Windows QTl cartographer 2.5. http://statgen.ncsu.edu/qtlcart/WQTLCart.htm. 
Churchiii GA, Doerge RW: Empirical threshold values for quantitative trait 
mapping. Genetics 1994 138:963-971. 
Phytozome. http/'/vmvw.phytozome.net. 



doi:10.1 186/1471-2229-14-52 

Cite this article as: Romeu ef a/.: Quantitative trait loci affecting 
reproductive phenology in peach. BMC Plant Biology 2014 14:52. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.blomedcentral.com/submit 



o 



BioMed Central 



